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ABSTRACT 

Electric currents j flow along the open magnetic field lines from the polar 
caps of neutron stars. Activity of a polar cap depends on the ratio a = j/cpcj, 
where pcj is the corotation charge density. The customary assumption a ~ 1 
is not supported by recent simulations of pulsar magnetospheres and we study 
^ . polar caps with arbitrary a. We argue that no significant activity is generated 

on field lines with < a < 1. Charges are extracted from the star and flow along 

O 

^ , is generated, leading to unsteady discharge on a scale-height smaller than 

the size of the polar cap. The discharge can power observed pulsars. Voltage 
^ . fluctuations in the discharge imply unsteady twisting of the open flux tube and 



O 



^ , generation of Alfven waves. These waves are ducted along the tube and converted 

CN ! to electromagnetic waves, providing a new mechanism for pulsar radiation. 

0^ ■ 

o ■ 

^ , Subject headings: plasmas — stars: magnetic fields, neutron 



Introduction 



X 

^ \ Corotation of a plasma magnetosphere is impossible beyond the light cylinder of a star, 

and magnetic field lines that extend beyond this cylinder are twisted, V x B 7^ 0. Thus, 
currents = (c/47r)Vx B are induced in the open magnetic flux tubes that connect the star 
(its "polar caps") to the light cylinder (Sturrock 1971). These currents are approximately 
force-free and flow along the magnetic field B. A basic question of pulsar theory is what 
voltage develops along the open tube to maintain these currents. It determines the dissipated 
power and creation that feeds the observed activity of pulsars. 

The customary pulsar model assumes that the electric current Jb extracted from the 
polar cap nearly matches cpcj, where pcj = ■ B/27rc is the corotation charge density 
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(Goldreich & Julian 1969). The deviation of current from cpcj was calculated as an eigen 
value of an electrostatic problem and found to be small (Arons & Scharlemann 1979). This is 
in conflict with recent global models of pulsar magnetospheres, which report {Jb — cpcj I ~ Jb 
(e.g. Contopoulos et al. 1999; Spitkovsky 2006; Timokhin 2006; see Arons 2008 for a review). 
A significant mismatch between Jb and cpcj can be expected on general grounds (Kennel et 
al. 1979). Currents Jb are determined by the magnetic- field twisting near the light cylinder, 
while pqj is a local quantity at the polar cap that is practically independent of Jb- 

The open tube is surrounded by the grounded closed magnetospher^ and may be 
thought of as a waveguide, filled with magnetized (ID) plasma. Compared to usual plasma- 
filled waveguides, it has two special features: (1) Current Jb is imposed on the tube. 
The twisted tube extends into the star, which is a good conductor and maintains Jb- 
sign(jB) = ±1 corresponds to ± charge flowing outward along the magnetic field lines. 
(2) Vacuum has effective charge density po = "~Pgj as Gauss law in the rotating frame reads 
V ■ E = 47r(p — Pgj)- The key dimensionless parameter (which can vary along and across 
the tube) is 



In this Letter, we discuss basic properties of the polar-cap accelerator with arbitrary a. 
First, we discuss what happens without creation: §[2] studies how the current is extracted 
from the polar cap of a radius and flows at small heights z <^ r^^ (this region is called 
"slab zone" below), and § [3] discusses how the flow extends to the region z > Vpc ("thin-tube 
zone" ) . We argue that the accelerator is inefficient if < a < 1 . 

The value of a depends largely on the angle x between f2 and B at the polar cap. For 
aligned dipole rotators {x = 0), magnetospheric models predict a < 1 everywhere on the 
polar cap and a < near its edge (see Fig. 4 and 5 in Timokhin 2006). For orthogonal 
rotator (x ^ 7r/2), |a;| ^ 1 throughout most of the polar cap. Generally, the polar cap has 
three regions where a;>l,0<a<l, and a < 0. We propose that pulsar activity originates 
in the polar-cap regions where < 1 (i.e. a > 1 or a < 0) as a high voltage is generated 
in these regions. 

We emphasize that the voltage is generated because V x B 7^ 0, not because pcj 7^ 0. 
The accelerator works as well if pcj = (i.e. if ± B at the polar cap), which corresponds 
to a — ±00. The accelerator height h ^ rpc is regulated by unsteady discharges (§ H]). 
§ [5] describes a mechanism for radio emission from unsteady discharges. 



^ The closed magnetosphere with Jb = is expected to have p ~ pG,i and « (e.g. Arons 1979). 
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2. Acceleration of a charge-separated flow at z <^r 



pc 



Suppose the current Jb is carried by charges e extracted from the polar cap; e < 
if Jb < 0. We assume that this charge-separated flow is freely supplied by the star with 
initial velocity v <^ c. Let us look for a steady state in the rotating frame of the star. 
The electrostatic fleld E is vertical near the surface (except at the edge of the polar cap). 
Since the displacement current c^^dE/dt = 0, the conduction current j satisfles j = = 
(c/Att) V X B. All particles are in the ground Landau state and flow along the magnetic 
fleld lines. Their motion is governed by the electric- fleld component E\\ (parallel to B), 



dt mc ' 

where m is the charge mass and p is its momentum in units of mc. Gravity is negligible here 
(in contrast to a similar problem for closed fleld lines, Beloborodov & Thompson 2007). 

Let ip be the angle between B and the vertical z-axis. Using V ■ E ^ dE^/dz for 
the electrostatic fleld near the surface {z <^ '"pc), one flnds from the Gauss law, dE\\/dz ~ 
47r cos ijj {p — pgj). The flow has charge density p = j/v where v = cp{l +p'^)~^^'^. In a steady 
state, one can rewrite Gauss equation using the convective derivative d/dt = cosip v d/dz, 



where t = corresponds to 2; = 0; -E||(0) ^ on the conducting surface. The parameter a{z) 
is not constant, because the magnetic fleld lines are curved (Scharlemann, Arons, & Fawley 
1978) and B is changed by the frame dragging effect (Muslimov & Tsygan 1992). At small 
z one can use the linear expansion a{z) ~ ao + ( z/R, where R is the star's radius. We 
will illustrate with ^ = — 1. Equations ([2D and ([3]), are solved together with the equation 
dz/dt = cp(l +p^)~^/^. The obtained steady-state models p{z) are summarized in Figure 1. 

One can understand the results by considering the simplifled model with a{z) = const, 
which describes a relativistic diode (e.g. Mestel et al. 1985). Let us also assume ip = 0. The 
flrst integral of equations ([2]) and is then given by 



where po = p{0) ^ 1. E\\ vanishes ii p = po oi p = 2a/ (1 — a^) + C(po)- The second root is 
positive if < a < 1; then p{z) oscillates between po and 



dp eE\ 




(3) 




2a 



(5) 
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a-' > 1. 
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Fig. 1. — Steady-state solutions for the charge flow extracted from the polar cap by the 
self-consistent E\\ in the slab zone z <^ rpc. The magnetic fleld lines are assumed to be 
perpendicular to the polar cap [tp = 0). Besides a{z), the problem has the dimensionless 
parameter ^ = {AtcJ e/mc^Y ^'^Vpc ~ {2u!bY^'^ ^ {R/ c)^^"^ where ujb = eB/mc. We assume 
in the flgure ^ = 10^ and Vpc/R = 10~^, which approximately corresponds to pulsars with 
P ~ 1 s and B ~ 10^'^ G at the polar cap. (In the special case |a| -C 1 one would need to 
reduce ^ by a factor ~ |a|, but we do not plot such cases here.) The oscillatory solutions 
(solid curves) will extend to the thin-tube zone z > rpc. The change in a{z) = ao + (z/R 
{( = —1 in this example) leads to ,20(2;) 7^ const and Pma.x{z) 7^ const; this effect is visible 
for curves with |a — 1| < \C\rpc/ R ~ 10~^ \C\P~^^'^. The monotonic solutions (dotted curves) 
illustrate that a high voltage is generated if a"^ < 1. They do not include creation 
that leads to intermittent screening of the accelerator (§ 4). Besides, the solutions neglect 
V_L ■ E± in the Gauss equation; this term would become comparable to j/c in an unscreened 
accelerator with the maximum p ~ ~ 10®. 
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The flow starts at E\\ = 0, with the minimum p = po ^ 1 and the maximum p/pcj ^ 1- 
Maximum -E^^x — Snjmce^^^l — Vl — a^) is reached where p/pcj passes through 1 {v = a). 
For < a ^ 1 and po <^ a, the solution p{t) (of eqs. [2] and H]) is a sinusoid p(t) = 
a (1 — coscUpt), where = inepcj/m. Using z{t) = J vdt one finds p{z). For < a <^ 1, 
p{z) is a cycloid with period zq = 2Trac/ujp. For 1 — a <^ 1, the period zq = 2'^/^(l — a)^^c/c<jp. 

By contrast, flows with a"^ < 1 (i.e. a > 1 or a < 0) quickly asymptote to 

In such flows, p does not approach pcj; instead it quickly saturates at j/c. Then edE\\/dz = 
47r(l — a~^)ej/c > 0, so eE\\ keeps growing oc z, and hence p grows oc z'^. 

The special flow with a = 1 was studied in detail (Michel 1974). The flow develops p > 1 
at 2; > c/ujp, and one finds -Ey = (87rmcj/e)^/^[l — 0{l/p)] = const and p{z) = \/2 ^z/r^c- 
This simple solution assumes a{z) = 1 = const. It was modified and extended to 2; > r^^ by 
Fawley et al. (1977); subsequent works studied the effects of slow variations in a near unity. 

3. Extension to the thin-tube zone z > and time dependence 

We still discuss here flows with no creation. The effects of are addressed in § |H 
Consider a solution from § |2]with a^^ > 1 and period Zq <^ r^^ (Fig. 1). The flow returns 
to its initial conditions at 2; = after passing distance Zq, i.e., effectively, the location of 
the polar-cap surface is shifted to Zn = uzq {n = 1,2,...). Therefore, the flow remains in 
the slab regime and its periodic acceleration/deceleration continues at z > Vpc with a small 
oscillating potential dropE e^/rrieC^ = (1 +Pmax)^^^ 1 = 2a^(l — a^)~^. The parameter 
a and $ gradually change on scale ~ radius r ^ zq. Thus, a modest voltage sustains the 
charge flow along the bundle of field lines with < a < 1. Energy is released with rate ~ /$ 
where I is the current along the bundle. It cannot feed pulsar activity because of small 

A time dependence does not qualitatively change the character of the flow with < a < 
1: p{z,t) remains oscillatory because E\\ always directs the flow toward v/c = a. Indeed, 
acceleration to v/c > a implies either p/pcj < 1 or j/js > 1- In both cases, E\\ is generated 
to decelerate the flow: p/pcj < 1 implies edE\\/dz < (Gauss law), and j/js > 1 implies 
edE\\/dt < (self-induction). This negative feedback leads to the oscillatory behavior. E.g. 



^ This expression is invalid if zq ~ '"pc i.e. if 1 — a ^ A = 2^/^(c/a;prpc) ^ 6 x 10 ^ P Bi^^^{m/me)^^'^. 
A formal maximum p = A^^ (with a fine-tuned to 1) is an overestimate as it neglects the change in a{z). 
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perturbations with wavevectors k || B are simple Langmuir oscillations. Oscillations of 
amplitude -Emax can nonlinearly interact with the plasma flow, heat it, and lose coherence. 
Turbulent pulsations of electric field can introduce anomalous resistivity and increase voltage. 
This, however, still does not lead to strong particle acceleration needed for pulsar activity. 

The basic difference between the regimes < 1 and > 1 is due to the interplay of 
self-induction and electrostatic effects. Consider first the case a = (zero current). Clearly, 
neutrality p = pqj should be maintained; the needed charges will be pulled out from the star 
and a static atmosphere can form with p ^ pqj. For a small a > 0, the charges will slowly 
drift outward and will be replenished by new charges extracted from the polar cap by a small 

This picture can remain qualitatively the same with increasing a until a approaches 1. 
When a > 1, the extracted charges with density p = pqj are unable to carry the required 
current Jb even if they move with v = c; then a large will be induced. In essence, it is 
the vacuum charge density po = — Pgj that helps avoid high voltages for < a < 1: then 
the extracted plasma both neutralizes the tube and carries the imposed current. Such an 
accelerator is inefficient because E\\ is dynamically screened on the small plasma scale 2:0- 
This general argument should hold in full 3D time- dependent models of the open tube. 

Hereafter, we focus on the regime a^^ < 1 (a < or a > 1), when strong particle 
acceleration is expected (Fig. 1). In absence of creation, the flow keeps accelerating until 
it approaches z ~ Vpc where eE\\ drops. Consider a simple case where a = const across the 
polar cap. If one imagines a global steady state with j = Jb, the electrostatic potential at 
z > would be (e.g. Ruderman & Sutherland 1975), 



where S± is the cross section of the open tube, a^^ < 1 implies e$ > 0, and the acceler- 
ated flow cannot overcome the potential barrier at the entrance to the thin-tube zone. It 
would have to decelerate and reverse at z ~ (in an attempted steady state, the deceler- 
ation would increase density, making the potential barrier higher). Therefore, j cannot flow 
steadily. When J/Jb drops below unity, the induction current {4:TT)~^dE\\/dt = Jb — j is gen- 
erated, and e£^|| grows until j recovers back to j = Jb- Thus, in absence of creation, the 
charge-separated flow from the polar cap is intermittently reversed by the electrostatic bar- 
rier at z ~ and then pushed through by the self-induction effect. The particle momentum 
p achieved in this unsteady regime is Pmax ~ (1 ~ -) (eq- E]). 




(7) 



- 7- 



4. Polar-cap accelerator with e discharges 



Accelerators with a < or a > 1 are efficient and can easily ignite an discharge that 
screens £"11 at a height h < r^c- The voltage developed in the unscreened "gap" < z < h is 



The gap with a > 1 was discussed by Sturrock (1971) and Kennel et al. (1979). It is 
different from the vacuum gap of Ruderman & Sutherland (1975), which formed because 
the star was unable to supply charges to keep p = pcj. The gap described by equation ([8]) 
forms because the imposed current Jb oversupplies charge density {p/pcj > 1 when a > 1) 
or supplies it with opposite sign (p/pcj < when a < 0). The electric field in the gap is 
much stronger than in Arons & Scharlemann (1979) model. The maximal available potential 
drop $max ~ BVt^R^ /c^ is achieved if /i ~ r^c- In observed pulsars h < r^c- 

The discharge develops as the accelerated primary particles (extracted from the polar 
cap) emit copious 7-rays that convert to outgoing e^. The Lorentz factor 7± ~ 10^ — 10^ is 
much smaller than the primary 71 ~ eE\\z/mc^. Therefore, the created — e charges (opposite 
in sign to jb) are easily reversed by E\\ and flow back to the polar cap. The reversal length 
in the unscreened E\\ is Zrcv ~ (7±/7i) z. 

The backflow would be avoided only if the rate of production aX z = h jumps from 
to a high value, launching a steady dense outflow that screens E\\ on a scale dz^cv < 4ev 
This, however, cannot happen for the following reason. The maximum charge density that 
may be created in Sz^cr by outflow without a backflow i^ p± ~ —^^M.scr3B/c, where 
A^scr is multiplicity (number per primary particle) of created in 5-2scr- The screening of 
£"11 requires — 47r5zscrP± ~ E\\, and 5zscr < 4ev leads to condition A^scr > l±liij'pc/ihY- 
It requires a huge and Cannot be satisfledH One concludes that a strong backflow is 

inevitable in pair-producing accelerators with h < r^c- 

A steady-state assumption for the discharge with a strong backflow leads to contradic- 
tions. First note that, in a steady state, the reversal of — e implies edp/dz > 0. It steepens 
the growth of eE\\ (z) at h ^ z ^ Vpc instead of screening E\\ . Thus, all charges — e created at 
z < Tpc reverse in the unscreened £11, creating a backflow of a large multiplicity A^back ^ 1- 



To avoid the backflow, the outflow must be gently polarized, with a typical reduction of 7_e by a 
factor ^ 2; otherwise the dispersion A^^^/l-e = would help form a strong backflow. 

Sharpest pair formation fronts occur in accelerators with 71 > 10^, where are created by curvature 
photons (e.g. Arons 1979). Other channels of pair production (significant at smaller 71) create e"^ with a 
smooth distribution over z and help backflow formation. 




h < r. 



pc- 
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The backflowing charges reach the polar cap and contribute a fraction / to the total current, 
leaving the fraction (1 — /) to be carried by charges extracted from the polar cap. In a 
steady state, the acceleration of extracted charges would be described by the same equations 
of § [2] if one substitutes j = (1 — /) and replaces pcj by pcj + fjB / c (the effective vacuum 
po now includes the backflow charge density —fjB/c). The new parameter aes is given by 
«efT = j/(cpGJ + fjs) = «(1 ~ /)/(! + ct/)- The presence of backflowing charges (/ 7^ 0) 
changes the condition for gap formation from a^^ < 1 to a~f^ < 1, which is equivalent to 

/</.4(i-i). m 

The charge density in the relativistic accelerator is Pgap ~ /js/c — (1 — f)jB/c — pcj = 
{l — 2f — a~^)jB/c. It must satisfy edE\\/dz = 47repgap > 0, which explains the condition ([9]). 
One then finds: (1) In a steady accelerator with a > 1, f cannot exceed ^. It permits 
a backflow multiplicity A^back = //(I ~ /) < 1- A steady discharge with a > 1 is 
inconsistent because it would create A^back ^ 1 that switches off the polar-cap accelerator. 
Similarly, a steady state is excluded for a < —1. (2) The accelerator with — 1 < a < 
could, in principle, have / — 1 and Alback ^ 1. A steady state is, however, implausible for 
another reason: a significant rate of ingoing pair creation can be expected when A^back ^ 1- 
It would lead to a runaway loop of creation in the accelerator, breaking the steady state. 
(3) A special case arises if / — > which requires A^back = A^back = (a — !)/(« + !)• Note 
that it cannot happen if — 1 < a < 1. For other a, Alback ~ ^ (unless a is fine-tuned to 
— 1). The accelerator should not permit a steady state with Alback = Alback ~ ^ because 
the discharge creates Alback ^ 1 even if f —> (details will be given elsewhere). 

We conclude that the discharge is unsteady, repeating on a timescale ~ /i/c. Fol- 
lowing a discharge, / exceeds f^, and voltage is screened. Voltage is re-generated when 
leave the polar-cap region and / drops below In contrast to steady models, E\\ is not 
described by an electrostatic potential $. Its time evolution obeys the Maxwell equation, 
dE\\/dt = 4:TT{jB — j), which expresses the self-induction effect. The current will fluctuate 
around Jb (Levinson et al. 2005); in addition, the excitation of Alfven waves implies fluc- 
tuations of Jb itself (§5). A time-averaged Alback ~ 1 can be expected in such discharges. 
Then the backflow energy flux is comparable to that of the primary particles. The resulting 
strong heating of the polar cap can conflict observations, posing a problem (cf. Arons 2008). 

5. Alfven- wave pulsar 

E 7^ in the open tube implies that the tube rotates relative to the polar cap (e.g. 
Ruderman & Sutherland 1975). The rotational velocity of magnetic field lines is given by 



-9- 



|v_l|/c = |Ex B|/S^ = E±/B where E <^ B and E± is the component perpendicular to 
B. v± is the drift velocity of plasma; it may also be associated with rotation of magnetic 
energy since cE x B/47r is the Poynting flux. 

Consider now a time-dependent state, where the charge density and E± fluctuate on a 
vertical scale h and timescale t ~ h/c. Such fluctuations will be accompanied by changes in 
v±, and the magnetic field lines will be twisted with amplitude 

5B^ 5vit 5vi_ SE± 

-B-^-h T ^ 

The corresponding fiuctuation in magnetic energy density is 5{B^/8n) = 5B • B/47r + 
{5 B)^/87r ~ {5B^Y /An. The time-dependent twisting generates Alfven waves that propagate 
along the open tube with nearly speed of light and carry luminosity 
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LA-^rl^c'-^--^rlA^E^Y. (11) 

For a gap of height h < Vpc with charge- density fluctuations dp, one can estimate 5E± ~ 
And ph. This gives 

PgapjB 

where Pgap — p — Pgj ~ ^/Anh?, Lg^p — is the time-averaged power dissipated in the 
gap, and / ~ S±jB- Lx should not exceed i^gap, since it is the gap accelerator that provides 
the energy for shaking the field lines. A related estimate of electro-magnetic dipole radiation 
from sparks is found in Fawley (1978). 

The generated Alfven waves are ducted along the field fines and eventually suffer Lan- 
dau damping (e.g. Lyubarsky 1996). However, part of La will convert to transverse elec- 
tromagnetic waves in the emission region, where the modes do not adiabatically track. Jil, 
Lyubarsky, & Melikidze (2004) show that such conversion generally occurs for waves with 
frequency below the plasma frequency. This condition is satisfied for Alfven waves and a 
significant fraction of luminosity La is expected to escape as radio emission. 

Usually considered mechanisms of pulsar radio emission are based on the two-stream 
instability in the outflow at large radii. The emerging radio spectrum is then determined 
by the plasma physics and possibly by curvature emission from plasma bunches created by 
the two-stream instability. By contrast, the emission fed by Alfven waves from the gap will 
have a spectrum controlled by the discharge behavior. One characteristic wavelength may be 
associated with the gap thickness /i, roughly corresponding to frequency ~ 10 MHz. However, 
fluctuations 5p will occur also on scales X < h, producing waves of higher frequencies. Such 
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fluctuations can create 5E± ~ 47r 5p A, and their contribution to Alfven luminosity is 



A short scale available in the discharge is the screening scale Zrev ^ h (§111). The electric 
field is screened by sudden formation of dense fronts of — e backfiow on scale A ~ /rev, pro- 
ducing large fiuctuations Sp/pcj 3> 1, which can generate waves of frequency c/Zrev ~ GHz. 
The spectrum of fiuctuations in the discharge extends from ~ c//i to c//rev Numerical 
simulations of the discharge can help find the spectrum of escaping waves and their beaming. 

This work was supported by NASA grant NNG-06-G107G. 
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